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Section  i.  INTRODUCTION 


This  report  is  concerned  with  the  problem  of  estimating  the  survival 
of  elements  of  a  penetrating  attack  force  using  a  variety  of  penetration 
aids.  This  problem  has  a  long  history,  having  been  investigated  by  a 
number  of  aircraft,  research,  and  government  organizations.  Most 
approaches  to  the  problem  are  based  on  the  fact  that  the  status  of  the 
penetrating  force  depends  on  a  time  sequence  of  events  with  the  property 
that  each  event  depends  upon  the  history  of  previous  events,  all  actions 
that  have  been  taken,  and  random  outcomes  associated  with  these  events 
and  actions. 


The  level  of  detail  in  the  modeling  usually  dictates  the  proper 
approach  to  the  problem.  Simple  formulations  may  be  handled  using 
elementary  methods  of  probability.  Complex  versions  of  the  problem 
require  either  Monte  Carlo  or  Markov  Chain  solutions.  Both  Monte  Carlo 
and  Markov  Chain  approaches  keep  track  of  changes  in  the  state  of  the 
penetration  process  as  discrete  time  events  and  random  outcomes  based 
on  these  events  and  the  current  state  of  the  process  occur.  The 
Markov  Chain  solution  actually  requires  discrete  time  updating  of  the 
probability  distribution  defined  on  the  state  space.  There  are  certain 
size  state  space  problems  for  which  the  Markov  Chain  approach  is  pre¬ 
ferred  over  the  Monte  Carlo  approach.  ^Generally  speaking,  if  the  size 
of  the  state  space  is  large,  Monte  Carlo  methods  are  preferred.  In 
penetration  studies  it  is  not  hard  to  come  up  with  problems  that  strain 
the  capacity  of  computers  and  the  budgets  of  those  performing  the  studied. 


This  report  offers  another  method  of  handling  penetration  problems. 
Penetration  problems  are  formulated  as  continuous  stochastic  processes 
using  concepts  from  the  area  of  pure  death  processes.  Whereas  time  is 
the  natural  parameter  of  interest  in  death  processes,  depth  of  penetra¬ 
tion  is  the  most  suitable  parameter  to  use  in  penetration  problems.  The 
event/probability  updating  concept  in  discrete  models  is  replaced  by  a 
penetration  interval/probability  updating  concept  in  continuous  models. 


In  this  modeling  the  solution  of  stochastic  differential  equations 
giving  the  probability  distribution  of  survivors  as  a  function  of 
penetration  depth  is  the  primary  tool  of  analysis.  In  some  cases  the 
stochastic  differential  equations  describing  the  evolution  of  the 
penetration  process  may  be  solved  in  closed  form.  Even  in  the  most 
complex  situations  numerical  procedures  for  solving  the  appropriate 
stochastic  differential  equations  are  conceptually  easy  to  implement; 
however,  computation  can  still  be  a  problem.  This  report  presents  a 
method  of  handling  the  complexities  of  synergistic  penetration  aid 
effects  while  maintaining  a  reasonable  level  of  computational  effort. 
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This  report  consists  of  four  sections  and  one  appendix  in  addition 
to  the  introduction.  Section  2  discusses  simple  discrete  penetration 
models  and  introduces  a  discrete  approximation  which  has  a  well-defined 
continuous  analogue  which  is  introduced  in  Section  3.  Numerical  calcu¬ 
lations  which  compare  the  discrete  models  with  the  discrete  approximation 
introduced  in  Section  2  are  found  in  Appendix  A. 

Section  3  discusses  the  basic  stochastic  processes  involved  and 
presents  general  model  structures.  Model  complexity  is  discussed  and 
various  methods  of  solving  the  stochastic  differential  equations  associated 
with  the  penetration  problem  are  given.  Generally  speaking,  one  should 
be  able  to  put  together  solutions  not  requiring  numerical  integration 
of  the  differential  equations.  The  general  model  for  handling  penetra- 
tors  of  different  types  is  presented,  and  statistical  independence  of 
the  penetrator  types  is  discussed. 

Section  4  is  devoted  to  the  calculation  of  the  kill  potential  that 
the  defense  can  bring  to  bear  on  a  penetration  force.  It  is  assumed 
that,  at  the  very  least,  the  attributes  of  a  single  site  against  a  single 
penetrator  are  known.  This  data  is  then  input  to  a  simple  calculation 
to  yield  cumulative  kill  potential  for  a  nonhomogeneous  defense.  The 
defense  is  also  allowed  to  restrict  firings  to  those  having  an  acceptable 
yield  in  terms  of  kill  probability.  Another  important  consideration 
in  defense  capability  is  the  allocation  of  defensive  firings  to  particu¬ 
lar  types  of  penetrators  and  the  fact  that  the  defense  may  make  errors 
in  the  classification  process.  If  data  of  a  single  site  as  a  function 
of  force  size  is  known,  this  data  may  be  input  to  the  model  directly. 

Section  5  is  devoted  to  the  effects  of  penetration  aids  and  tactics. 

In  this  section  cumulative  kill  potential  calculations  are  modified  to  1 

account  for  jinking  and/or  maneuvering,  ECM,  defense  suppression,  self- 
defense,  and  standoff  tactics.  Since  the  model  is  structured  to  handle 
nonhomogeneous  penetration  forces,  it  is  a  simple  matter  to  include 
decoys  in  the  model.  The  analysis  of  synergistic  effects  in  penetration 
aid  studies  can  lead  to  extremely  large  state  space  problems.  An  iterative 
procedure  for  including  synergistic  effects  while  maintaining  a  reasonable 
size  state  space  is  discussed. 
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Section  2.  SIMPLE  DISCRETE  PENETRATION  MODELING 


In  this  section  we  want  to  discuss  simple  discrete  penetration 
models.  We  also  want  to  consider  a  discrete  approximation  to  these 
discrete  models  which  has  a  continuous  analogue.  This  continuous  model 
is  derived  from  basic  principles  in  Section  3.  We  assume  that,  with 
each  defense  site,  there  is  a  maximum  number,  M,  of  independent  defense/ 
penetrator  engagement  opportunities  which  the  defense  can  mount  against 
a  penetrating  force.  The  parameter  M  reflects  all  the  limitations  that 
are  imposed  on  the  defense  site  by  recycle  times,  penetration  time  in 
envelope,  switching  times,  ammunition  limitations,  and  reload  times. 

For  now,  we  assume  each  engagement  has  a  probability  of  kill  of  Pk. 

An  engagement  may  imply  the  launch  of  a  salvo  of  some  number  of 
missiles  or  the  use  of  some  nominal  number  of  AAA  rounds  and  then  Pk  is 
well  defined  with  respect  to  the  definition  of  engagement  for  a  given 
defense  site.  The  most  common  definition  of  engagement  will  be  the 
launch  of  a  single  missile  or  the  firing  of  sufficient  AAA  rounds  to 
yield  a  prescribed  kill  probability,  Pfc.  In  the  latter  instance,  a 
common  situation  is  M  =  1  and  then  Pfc  is  a  consequence  of  continuous 
fire . 


2.1  PENETRATOR  GROUP  SIZE  N  >  M 

When  N,  the  number  in  the  penetrating  group  is  at  least  as  large 
as  M,  one  is  led  to  a  simple  solution  for  the  number  of  penetrators 
killed.  In  this  case,  the  defense  might  as  well  utilize  all  M  oppor¬ 
tunities.  The  probability  that  i  out  of  the  M  penetrators  engaged  are 
killed  is 

P(i)  =  (i)  Pk(l  '  Pk)M_i  i  =  0,1,*  •  •  ,  M  (2.1) 

with  expectation 

E  =  HPk  (2.2) 

Frequently  one  wants  a  definition  of  defense  capability  that  is 
independent  of  the  size  of  the  penetration  force.  A  common  measure  of 
defense  capability  is  the  product  of  engagement  opportunities  and  kill 
probability  P^  defined  to  be  the  kill  potential, 
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(2.3) 


In  this  particular  instance  when  N  £  M,  E  =  H.  The  kill  potential  of 
a  site  is  a  particularly  simple  measure  of  defense  site  effectiveness 
and  it  arises  naturally  in  certain  discrete  approximations  to  the 
binomial  distribution  Equation  2.1. 

In  fact,  if  M  -*■  00  and  -*■  0  in  such  a  way  that  MPfc  *  H,  then 
Equation  2.1  may  be  approximated  by  the  Poisson  distribution  giving 


P(i)  =  e 


-H  H_ 
i! 


i  =  0,1, 


(2.4) 


The  implication  is  that  N  -*■  00  since  N  k  M  and  this  approximation  will  be 
valid  for  all  i  only  in  those  situations  in  which  N  and  M  are  both  large 
and  P^  is  small. 

2.2  PENETRATOR  GROUP  SIZE  N  <  M 

When  N  is  less  than  M,  we  want  to  consider  two  possible  defensive 
actions.  One  possibility  is  to  allocate  multiple  simultaneous  engage¬ 
ments  (MSE)  to  a  single  penetrator.  Salvo  and  ripple  firing  are 
examples  of  allocating  more  than  one  engagement  at  a  time  to  a  single 
penetrator.  The  second  possibility  is  to  use  the  M  encounters  in  such 
a  way  that  an  encounter  is  always  allocated  to  a  penetrator  alive,  not 
currently  under  attack.  Shoot-look-shoot  (SLS)  firing  doctrines  fall 
into  this  category. 

Suppose,  for  example,  N  =  1,  H  =  2  and  there  is  sufficient  time  to 
employ  an  SLS  firing  doctrine.  The  probability  of  kill  of  the  penetrator 
is 

P  =  P,  +  (1  -  P.)  P. 


and  the  expected  number  of  kills  is 


E  =  P  =  2P.  -  P, 
k  k 


\ 


•JL 

.  J 
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Since  the  kill  potential  is 


H  "  2Pk  • 


E  <  H 


In  an  SLS  firing  doctrine  it  is  possible  to  kill  the  penetrator  with  a 
single  engagement.  If  an  engagement  corresponds  to  the  use  of  a  weapon, 
the  expected  weapons  usage  is 


Ew=lPk+  2(1  -  Pk)  =  2  -  Pk 


For  an  MSE  firing  doctrine  we  have  the  corresponding  results  desig¬ 
nated  by  an  asterisk, 


p*  - 1  -  a  -  V2  - 2  pk  -  pk  ‘ F 


E*  =  p*  =  P  =  E 


with 


E*  <  H 


In  this  case,  however,  we  always  use  two  weapons,  so 


E*  =  2  >  2  -  P,  =  E 
w  k  w 


More  generally,  for  an  SLS  firing  doctrine  the  probability  that  i 
penetrators  are  killed  is 


P(i)  - 1")  pjci  -  pk>’ 


i  =  0,1,* • • jN-1 


-  E  (!)  - 


.  .  •  TJ,  ■*/  .... 


E 


-<  E  i(!)  pk<1  -  v"'1  -  k 


So  that  an  SLS  firing  doctrine  leads  to 


whenever  N  <  M.  The  expected  weapon  usage  is 


M-N-l 


E„-  E  (N  +r)(N+p-1)p“(l- 


“  E  (NTl)  -  vr 


r=M-N 


<  N  +  E  r  (N+p-1)  p“(l  -  Pk)r 
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In  this  case, 


SMBkSMNMt 


Clearly  E  <  M,  so 
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E  <  min 
w  - 


f>L  M 

K ' 


For  an  MSE  firing  doctrine  we  assume  M/N  is  an  integer  and  this 
number  of  engagements  is  placed  on  each  of  the  N  penetrators. 
Immediately 

E*  =  M  >  E 
w  w 

with  strict  inequality  whenever 

N  <  M  P. 

k 

The  probability  that  i  penetrators  are  killed  is 


p*(i)  =  (j)(i  -  (i  -  a  -  v 


.M/N (N-i) 


with  the  expected  number  killed  given  by 


=  n(i  -  (i  -  Pk)M/N) 


One  may  show  that  E*  is  an  increasing  function  of  N,  so 


<  M(l  -  (1  -  Pk)M/N)  =  M 


Pk  =  H 


It  is  fairly  easy  to  see  that  the  expected  kills  in  an  SLS  firing 
doctrine  will  exceed  those  obtainable  from  an  MSE  firing  doctrine. 

Suppose  the  MSE  doctrine  allocates  n^  engagements  to  the  ith  penetrator 
and  an  SLS  doctrine  is  used  on  each  penetrator,  stopping  if  the  penetrator 
is  killed  or  if  n^  engagements  are  utilized.  For  the  ith  penetrator  the 
probability  of  kill  is 


p'jffiwii  nun 


.....  . 


.•**»*-  .  .  .  >:  / 


m 
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P  =  Pk  +  \  Pk  +  “•  +  ql 


(■  -  .;■) 


1  -  <U 


1  -  <u 


which  is  the  probability  of  kill  of  the  salvo  of  size  n^  in  the  MSE  firing 
doctrine.  In  this  modified  SLS  doctrine,  providing  P^  >  0,  there  is 
positive  probability  that  not  all  n^  engagements  will  be  used  for  all  i. 

It  is  also  clear  that,  providing  P^  <  1,  there  is  positive  probability 
that  at  least  one  penetrator  survives.  If  6  is  the  probability  that 
not  all  engagements  are  used  and  not  all  penetrators  are  killed,  then 


E  >  E*  +  6  P,  >  E* 
k 


since  in  an  SLS  firing  doctrine  the  remaining  engagements  could  be 
allocated  to  any  survivor. 


2.3  A  POISSON  APPROXIMATION 


It  was  observed  in  2.1  that  under  certain  instances  the  Poisson 
distribution  may  be  used  as  an  approximation  to  the  binomial.  In 
general,  M  may  be  large,  but  P^  may  or  may  not  be  large.  Generally, 


(?)  pt(1  -  pk>M'i  *  e 


-MPk  (M  Pk>: 


for  M  (depending  on  i)  large  enough,  regardless  of  P^.  If,  moreover, 
N  <<  M,  then  we  are  concerned  with  the  accuracy  of  the  above  approxi¬ 
mation  for  i  =  0,1,* ■•jN. 


These  comments  suggest  the  following  approximation  for  the  binomial 
probabilities  under  consideration. 


HMH '  w  i  *  "PTy 
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P  (N) 


where 


H  =  M  P, 
k 

We  shall  not  specify  the  relationship  of  M  to  N,  but  one  feels  that  for 
large  M  the  approximation  will  be  adequate. 

In  the  next  section  we  shall  see  that  reasonable  assumptions  concern¬ 
ing  the  attrition  of  penetrators  under  attack  leads  to  a  stochastic 
process  in  which  the  number  of  penetrators  killed  1  as  a  truncated 
Poisson  distribution.  We  shall  use  cumulative  kill  potential  of  defense 
installations  the  penetration  force  encounters  as  the  parameter  in  this 
distribution.  Moreover,  it  is  our  intent  to  modify  kill  potential  to 
reflect  both  defense  and  offense  strategies  in  achieving  their  individ¬ 
ual  goals. 


In  this  section  we  have  seen  that  the  truncated  Poisson  can  be  good 
approximation  to  the  binomial  when  SLS  is  the  firing  doctrine  and  M  is 
large.  Appendix  A  provides  a  numerical  comparison  of  the  expectation 
of  the  truncated  Poisson  distribution  with  parameter 


I 


H  =  M  P, 
k 

with  the  binomials  associated  with  both  SLS 
when  M  is  not  necessarily  large. 


and  MSE  firing  doctrines 


V 


11 


NWC  TP  5841 


Section  3.  BASIC  CONTINUOUS  STOCHASTIC  PENETRATION  MODELS 


In  this  section  we  consider  a  number  of  penetration  models  that  may 
be  solved  in  closed  form  by  solving  an  appropriate  system  of  differential 
equations.1  In  Sections  4  and  5  we  modify  these  models  to  account  for 
defensive  allocation  of  effort  and  the  effects  of  penetration  aids  and 
tactics.  Selecting  an  appropriate  stochastic  model  is  not  routine.  The 
selection  depends  upon  an  assessment  of  how  the  defender  will  utilize 
his  resources  and  manage  his  defense.  The  basic  question  is  whether  a 
defender  utilizes  as  much  of  his  resources  as  possible  or  whether  he 
paces  his  expenditure  of  ammunition  so  that  he  will  be  in  a  position  to 
defend  against  following  attacks. 

3.1  A  BASIC  POISSON  ATTRITION  MODEL 

We  assume  that  N  identical  penetrators  begin  a  penetration  of  enemy 
territory  at  x  =  0.  In  penetrating  from  0  to  x  the  penetrating  force 
encounters  a  number  of  defensive  sites  of  the  same  type.  We  assume  there 
exists  a  function  A(x)  having  the  property  that,  no  matter  how  many 
penetrators  have  been  killed  over  (0,x)  and  providing  at  least  one 
penetrator  is  alive,  then  the  probability  that  a  kill  occurs  over 
(x,  x  +  Ax)  is  A(x)Ax  +  0(Ax)  where 


Lim  0(Ax)  =  _ 
Ax-K)  Ax  * 


Moreover,  the  probability  of  more  than  one  kill  over  (x,  x  +  Ax)  is  0(Ax) . 
If  N  penetrators  have  been  killed  over  (0,x)  the  probability  of  an 
additional  kill  over  (x,  x  +  Ax)  is  0. 

These  considerations  lead  to  a  system  of  stochastic  difference  equa¬ 
tions  for  Pi(x) ,  the  probability  that  i  penetrators  are  killed  over 
(0 ,x) .  We  have 

Pq(x+Ax)  =  Pq(x) (l-A(x)Ax)  +  0(Ax) 

P_^(x+Ax)  =  P^(x) (l-A(x) Ax)  +  P^  ^(x)A(x)Ax  +  0(Ax) 

i  =  1 ,2 ,  • • • ,N-1 

PN(x+Ax)  =  PN(x)  +  PN_1(x) A(x)Ax  +  0(Ax) 

1  W.  Feller.  An  Introduction  to  Probability  Theory  and  Its 
Applications.  New  York,  John  Wiley  &  Sons,  Inc.,  1967.  Chapter  17. 
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Dividing  each  equation  by  Ax  and  taking  the  limit  as  Ax  -*■  0  yields 
the  system  of  stochastic  differential  equations, 


Pq(x)  =  -A(x)  PQ(x) 


P^(x)  =  -A(x)  P.(x)  +  A(x)  Pi_1(x)  i  =  1 ,2 , • • • ,N-1 


P^(x)  =  A(x)  PN_x(x) 


If  Pq(0)  =  1>  the  solution  of  this  system  of  differential  equations  is 


PjU)  -  e~H(x) 


i  =  0,1. •••.N-l 


where 


N-l 

VX)  =  1  -  E  Pi<x> 

i=0 


■/' 


H(x)  =J  A(y)dy 
0 


We  define  A(x)  to  be  the  kill  intensity  of  the  defense  at  x.  Based 
on  the  results  of  Section  2  and  Appendix  A,  we  assume  H(x)  gives  the 
cumulative  kill  potential  that  will  be  brought  to  bear  on  the  penetrating 
force  ever  (0,x).  We  note  that  if  Q^(x)  is  defined  to  be  the  probability 
that  i  penetrators  are  alive  at  x,  then 


Qi(x)  =  PN_i(x)  i  =  0 ,1, • • • ,N 

so 

N 

Qn(x)  =  1  -  ^  Qi(x) 
i=l 
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where 


Q±(X) 


=  iH00iN_i  -H(x) 
(N-i) 


i  =  1,2,  •  •  ■  ,N 


Methods  of  estimating  H(x)  as  a  function  of  defensive  resource 
management  policies  and  penetration  force  size  will  be  given  in  Section  4. 
Degradations  to  H(x)  due  to  penetration  aids  are  presented  in  Section  5. 
For  now  we  turn  our  attention  to  a  model  in  which  the  kill  intensity  is 
allowed  to  be  a  function  of  penetrator  force  size. 


3.2  A  MODEL  WHICH  ASSUMES  ATTRITION  IS 
PROPORTIONAL  TO  FORCE  SIZE 

One  reason  for  allowing  kill  intensity  to  depend  on  force  size  is 
the  simple  fact  that  a  reasonable  defensive  firing  doctrine  would  allocate 
more  kill  intensity  to  larger  forces.  After  all,  the  defense  must  balance 
available  resources  to  defend  over  extended  periods  of  time.  The  assump¬ 
tion  that  defensive  resource  allocation  will  be  related  to  the  size  of 
the  penetration  force  is  not  unreasonable. 

In  this  case  it  is  a  little  easier  to  set  up  the  differential 
equations  in  terms  of  the  number  surviving  over  (0,x).  If  i  penetrators 
survive  over  (0,x),  we  let  the  kill  intensity  at  x  be 


A^x)  =  i  <}>(x) 


where 


<Kx)  =  A(x)/N* 

where  N*  >  N  relates  to  the  defense  policy  of  allocating  kill  intensity 
to  various  sizes  penetration  forces.  This  guarantees  that 

A^Cx)  <  A(x)  i  =  0 ,1,2 ,  • •  •  ,N 


and  assumes  that  kill  intensity  is  allocated  proportional  to  the  size  of 
the  penetration  force.  By  setting  N*  =  N,  the  kill  intensity  when  none 
are  killed  is  the  same  as  in  3.1.  The  differential  equations  are 
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Qq(x)  =  4>(x)  Qi(x) 


Q^(x)  =  -i  <f>(x)  Q^x)  +  (i  +  1)  4>(x)  Qi+1  (x) 


i  =  1 ,2 , • • • ,N-1 


QjJj(x)  =  N  4>(x)  Qn(x) 


with  solution 


Q±(x) 


-(>- 


-$(x) \N  *  -i$(x) 

e  I  e 


where 


$(x)  =  /  <j>(x)dx  = 


X(x)/N*  =  H(x) /N* 


The  probability  distribution  of  the  number  killed  is  given  by 


Pi(x)  =  QN-i  (x) 


3.3  MIXTURES  OF  MODELS  IN  3.1  AND  3.2 

In  the  initial  phase  of  a  penetration  it  may  be  appropriate  to 
allocate  kill  intensity  proportional  to  the  force  size.  In  later 
phases,  however,  the  defense  may  be  willing  to  allocate  considerable 
kill  potential  to  even  small  forces.  This  suggests  the  possibility 
of  assuming  that  different  models  hold  over  different  subintervals  of 
(0 ,x) .  Suppose,  for  example,  the  model  of  (3.2)  holds  over  (0,X]_)  and 
the  model  of  (3.1)  holds  over  (x^,x).  In  this  case, 


-  (iX1  - 


-$(x)\N-i  -i$(x) 

e  I  e 


0  <  x  <  Xi 
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and 

N 

Qi(x)  =  £  Qj(xi)  Qi|  j(xi» x) 

j=i 

where 


Qi|j(xi' 


x) 


[H(xlf  x)]j_i 
_______ 


-H(x| ,x) 
e 


j 

Q0| j(xl»  x)  "  1  -  2  Qi|j(xl’  x) 
i=l 

and 


H(x1,  x) 


dy 


More  generally,  suppose  one  has  obtained  Qi(x)  over  the  subintervals 
(0,  ,  (x^,  X£) ,  (x^-l,  xk)  and  now  wants  to  obtain  Qi(x)  for 

x  >  xk.  In  this  case, 


Qx(x) 


Qj(xk>  ‘hi/v  x) 


X,  <  X 

k 


where  Qj (xk)  is  known  and  Qi|j(xk,  x)  depends  upon  the  model  assumed  to 
hold  over  x  >  xk.  If,  for  x  >  xk,  the  model  of  3.1  holds 


x» 


[H(x,  ,  x)]'’"1 

( j  -  i)  7 


-H(xk,x) 

e 


i  =  1,2,' 
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Q0|j(xk’  x)  =  1  "  2  Qi|j(xk’  x) 


H(xk,  x)  =  /  A(y)  dy 


If,  for  x  >  xk,  the  model  of  3.2  holds 


Qi|j(xk’  x)  =  (i)v1  -  e 


-<J)(xk,x)\^  1  -l$(xk,x) 


i  =  0,1,  • ■ •  ,j 


$(xk’  x)  =  /  X(y)  dy 

k+1  J 

xk 


where  Nk+^  >  j  is  a  defense  scaling  parameter  for  (xk,  x)  when  j  are  alive 
at  xk. 

3.4  A  COMBINED  MODEL 

Another  method  of  combining  the  two  models  assumes  that  the  kill 
intensity  at  x  when  i  penetrators  are  alive  is 


A^x)  =  to  ( x)  +  i<Kx) 


i  =  1 ,2 , • • •  ,N 


Aq(x)  =  0 


w(x)  +  N*tj)(x)  =  A(x) 


H«r 


One  may  proceed  solving  the  linear  differential  equations  in  turn  to 
obtain  P2»P3>  etc.  These  expressions  get  quite  complex  very  quickly 
and  it  may  actually  be  better  to  solve  these  equations  numerically. 
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3.4.1  A  Special  Case 


In  the  special  case  that  co(x)  =  U)  and  <j>(x)  =  independent  of  x 
the  system  of  differential  equations  above  may  be  solved  in  closed  form. 
One  may  verify  that  the  solution  is  given  by 


P0(x)  -  e-(“*tN*x) 


P.(x)  -  e-C<®K»-i)«(1^-*«)  fj 


i  =  1,2,- • • ,N-1 


PN(x)  =  1  ‘  5Z  Pi(x) 


This  closed  form  solution,  together  with  the  idea  of  combining  models  that 
are  different  for  different  subintervals  of  (0,x),  provides  a  reasonable 
method  of  solving  this  problem  in  general.  We  merely  approximate  to(x) 
and  <)>(x)  by  piecewise  constant  functions  with 


oo(x)  =  u)k 


*k-i 5  x  -  *k 


k  =  1 ,2 , • • •  ,m 


<Kx)  =  4>k 


Vi  -  *  -  xi 


k  =  1,2 , • • •  ,m 


where  Xq  =  0. 


Designating  the  above  probabilities  by  Pj(x;  w,  ,  N) ,  we  have 


Vx)  =  2  pj<xk-i)pi-j(x-xk-i;  v  v  N-j)  vi  - x  -  xk 

j=0 


for  k  =  l,2,*-«,m  v;here  Pq(0)  =  1  and 


P.(x)  =  P.(x;  U)n  ,  ,  N) 


0  <  x  <  x. 


*  .  pj WPR 
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These  formulae  are  easily  modified  to  account  for  reliability.  Suppose 
each  penetrator  is  subject  to  a  constant  failure  rate  8  =  yv,  where  v 
is  the  penetration  speed.  All  one  need  do  is  replace  <f>k  by  <}>k+y  in  the 
above  formulae.  If  the  penetration  force  travels  a  distance  a  before 
beginning  the  penetration  at  x  =  0  then 


i 


II 

* 

'w' 

P-4 

E  Vvi)pi- 

V  ^k+^‘  N-^ 

(_i. 

ii 

o 

xk_i  <  X  <  xk 

e 

CM 

where 

P.(0)  = 

P.(a;  0,  y,  N) 

i  =  0,1,2, ••• ,N 

= 

GM1--**)1  • 

-(N-i)ya 

Another  method  of  obtaining  an  approximate  solution  is  to  obtain 
the  average  of  uj  and  <p  over  (0,x). 


0)  = 


/ 


x 


w(y)dy/x  =  ft(x)/x 


and 

<f>  =  $(x)/x 

giving 


P0(x)  -  e-(“ <*)«♦<*» 

i i_1 

P  =  e-(^(x)+(N-i)$(x))|1_e-$(x))  J-J  (a(x)+(N-k)<Kx)) 


MMHi mmmmmmmh 
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N-l 

PN(x)  =  1  "  ]C  Pi(x) 
i=0 


3.4.2  The  Expected  Value  of  Number  Killed 
Realizing  that 


N 


dE(x) 

dx 


iP^(x) 


i=l 


where  E(x)  is  the  expected  number  of  penetrators  killed  over  (0,x),  we 
obtain,  after  some  algebraic  manipulation  of  the  differential  equations 
in  this  section. 


d^X-’)-  +  <})(x)E(x)  =  0}(x)(l-PN(x))  +  N<Mx) 


with  solution 


E(x)  =  e 


-<Kx) 


x 

ai(y)(l-PN(y))e$(y) 


dy 


+  N[j..e-*<x>| 


In  most  instances  Pjj(y)  -  ^  ^or  s0 


E(x)  *  e  $(x)  J  u(y)e$^y)dy  +  N[l-e  $(x)) 

0 

Thus,  if  one  is  only  interested  in  the  expected  number  killed,  a  reason¬ 
able  estimate  of  this  expectation  is  available. 
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3.5  NUMERICAL  SOLUTION  OF  DIFFERENTIAL-DIFFERENCE 
EQUATIONS 

We  assume  that  the  kill  intensity  at  x  given  that  i  penetrators 
are  alive  is  A^(x).  We  have  considered  the  case  where 

Xi(x)  =  oj(x)  +  i(<j>(x)+y) 

and  1J  is  a  reliability  factor.  We  now  allow  A^(x)  to  be  quite  general. 
The  basic  difference  equations  are 

Q^(x+Ax)  =  (l-Xi(x)Ax)Qi(x)  +  Ai+1(x)AxQi+1(x) 

i  =  0 ,1 , • • • ,N 


where  Aq(x)  =  0  and  Qjjq.^(x)  =  0  for  all  x  and 
Ql(0)  -  (“)e-^1(l-e->,a)N'1 


with  the  understanding  that 

v°>  - 1 

if  y  =  0  or  a  =  0. 

For  a  fixed  Ax  the  computation  is  straightforward.  It  is  clear 
that  both  accuracy  and  cost  will  depend  on  Ax  and  experience  with  the 
solution  should  dictate  a  choice  of  Ax  which  balances  accuracy  with  cost. 


3.6  PENETRATORS  OF  DIFFERENT  TYPES 

If  the  penetration  force  consists  of  J  different  types  of  penetrators, 
then  the  stochastic  difference  equations  are 


Q.  (x+Ax) 

,XJ 


J 

2J 

j=l 


J 

+  2  y*;1! 

j-i 


’Vl’' 


,iJ)AxQi 


1’ 


»i 


j+1* 


»i 


(x) 
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i 

1 


where  Xj  (x;i^,  •  •  •  ,ij)  is  the  kill  intensity  against  the  jth  type  of 
penetrator  at  x  when  ij,j=l,***,J  penetrators  of  type  j  are  alive. 

The  differential  equations  become 


J 


Xj(x;i^,i2»‘ 

j=l 


»iJ)Q1  i 
J  ^1 > 12  * 


»i 


J 


(x) 


J 

+  X^,  (x }  i-^ » i2  >  “ 

j=l 


j+1* 


’^■J^i  i 
J  ^i  ’  2  * 


If 


Xj  (x;i1 , 


•S' 


,ij)  =  (x) 

.1 


independent  of  i]_,  •  •  •  ,ij_l,ij+i>  *  *  *  »ij  f°r  all  j  »  then 


1 


I 

j# 


1 


«i- 


,i 


(x)  =  n  q.  (x) 


j-1 


where  Qj^(x)  satisfies 


Q|  (x)  -  -Xi  (x^  (x) 

j  j  j 


+  X  Q,  (x) 
±i+l  j+1 
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of  Q 


To  verify  that  this  is  the  solution  we  have  from  the  assumed  form 


ili 


.ij 


(x) 


J  J 


)i1,---,i,(x>  ‘  Z  77  <*K.oo 

1  J  i-l  "  J 


J  J 


-  -  Z  n  «±  <*> 

« _  i  •  in 


j=l 


A±  (x)Q.  (x) 
3  j 


J  J 


+  Z  77  Oi  <*> 

,  ,  m 


3=1 


A  (x)Q  (x) 
j+1  3+1 


“  Y]  X.  (x)Q 


(x) 


+  A  (x)Q 

j=1  Vi  ii,**’,ij+i,**',ij 


(x) 


This  is  an  important  result,  for  it  says  that  survival  of  different 
types  of  penetrators  in  a  penetration  group  may  be  treated  separately 
providing  the  attrition  rate  of  a  given  type  does  not  depend  on  the 
number  of  penetrators  of  other  types  that  are  alive.  As  one  might 
expect,  this  is  rarely  the  case  when  penetration  aids  and  uncertainties 
in  the  classification  of  target  types  are  incorporated  in  the  modeling. 

We  may  certainly  extend  the  numerical  solution  suggested  in  3.5  to 
the  case  where  there  are  a  number  of  penetrator  types  and  independence 
does  not  hold.  Our  concern  is  over  sheer  magnitude  of  the  number  of 
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f  "Vfc?.,-'.  w  ■ 


l 


i\ 


h 

H 


| 

8 


calculations  required.  Suppose,  for  example,  we  have  a  penetration 
force  of  five  types  with  nine  penetrators  per  type.  In  the  case  of 
independence,  50  probabilities  must  be  obtained  each  time  the  joint 
probability  distribution  is  determined  at  x  +  Ax  from  the  results  at  x. 
For  non-independence,  the  number  50  grows  to  100,000  indicating  the 
computational  difficulties  encountered  when  independence  does  not  hold. 


By  combining  expected  value  concepts  with  the  stochastic  differential 
equations  that  arise  when  independence  does  not  hold,  one  may  solve  a 
sequence  of  independent  problems  that  approximates  the  dependent  Ccise. 

We  address  this  problem  in  Sections  4  and  5. 
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Section  4.  CUMULATIVE  KILL  POTENTIAL  FROM  THE 
DEFENSE  POINT  OF  VIEW 

In  this  section  we  want  to  develop  some  simple  models  for  estimating 
cumulative  kill  potential,  H(x) ,  where  H’(x)  =  A(x)  is  a  basic  input  to 
the  models  discussed  in  Section  3.  An  input  to  this  calculation  is  the 
kill  potential  of  an  individual  site  and  considerations  of  how  that  site 
will  utilize  its  available  kill  potential.  Cumulative  kill  potential 
calculations  are  extended  to  include  different  types  of  defense  units 
and  the  problem  of  allocation  of  kill  potential  to  penetrators  of  dif¬ 
ferent  types  is  also  discussed.  These  considerations  generally  lead  to 
penetration  models  which  may  be  solved  in  closed  form  using  the  models 
of  Section  3. 

There  are  aspects  of  the  defense  problem  that  lead  to  penetration 
models  that  must  be  solved  numerically  or  involve  dependence  consider¬ 
ations.  One  problem  of  special  importance  is  the  effect  of  mistakes  in 
classification  of  penetrator  types.  Modeling  considerations  that  allow 
one  to  simplify  this  problem  are  discussed. 

4.1  KILL  POTENTIAL  OF  A  SINGLE  SITE 

We  consider  a  group  of  penetrators  at  altitude  h  and  offset  from 
the  defense  site  a  distance  y.  The  defense  can  get  off  M  shots  with 
Pk's  as  indicated  in  Figure  1. 


FIGURE  1. 
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Actually  M  and  the  Pk's  depend  on  h,  y,  and  v,  the  speed  of  the 
penetrators.  As  a  consequence,  the  kill  potential  of  the  site  is 


M(h,v,y) 

Kp(h,v,y)  =  ^  Pk(h,v,y,m) 


where  M(h,v,y)  is  the  number  of  shots  as  a  function  of  h,  v,  and  y  and 
Pk(h,v,y,m)  is  the  kill  probability  of  the  mth  shot  as  a  function  of 
h,  v,  and  y.  In  general,  y  is  considered  uniformly  distributed  over 
(0,  R(h))  where  R(h)  is  the  effective  range  of  the  site  as  a  function 
of  h.  Consequently,  we  set  the  kill  potential  of  the  site  equal  to 


y.R(h) 

vh>v)  ■  rW  J  yh'v'y)dy 


Another  consideration  of  importance  to  the  defense  is  expected  weapons 
usage 


W(h,v)  =  J  M(h ,v,y) dy  =  M(h,v) 


The  average  probability  of  kill  at  offset  y  is 
M(h,v,y) 

pk(h,v,y)  =  Pk(h,v,y,m)/M(h,v,y) 


yh’v>  =  roo 


/.R(h) 

J  M(h,v, 


y)pk(h,v,y)dy 
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and  we  define  Pk(h,v)  so  that 

K  (h,v)  =  M(h,v)Pk(h,v)  =  W(h,v)Pk(h,v) 

4.1.1  The  Effect  of  Salvo  Firing 

Should  the  defense  elect  to  fire  in  salvos  with  r  missiles  per 
salvo  then 


M(h,v,y) 

Pk(h,v,y;r)  =  ^  (1-(1-Pk(h,v,y,m))r)/M(h,v,y) 

m=l 

and 


>.R(h) 

K  (h,v;r)  -  I  M(h,v,y)Pk(h,v,y ;r)dy 

0 


and  expected  weapons  usage  is 
W(h,v;r)  =  rM(h,v) 


Again  wa  define  Pk(h,v;r)  so 

Pk(h,v;r) 

Kp(h,v,r)  =  M(h,v)Pk(h,v;r)  =  W(h,v;r)  - - - 

4.1.2  Selective  Salvo  Firing 

We  consider  the  permutation  ii»i2»  *  *  *  ’^(h.v.y)  of  inte8ers 
1,2, • • • ,M(h,v,y)  with  the  property  that 


Pk(h,v,y,i1)  >  Pk(h,v,y,i2),---,  *  Pk(h,v,y,iM(h>v>y)) 
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Next  we  define 

Ifi  =  Uj  :  j  =  1,2, 

A 

so  that  Ijij  indexes  the  integers  associated  with  the  M  largest  Pk's, 
where  M(h,v,y)  <  M(h,v,y) 

%  "  UJ  ;  PfcCh.v.y.ij)  >  Pt) 

where  this  set  may  be  empty.  Thus  Ip,  identifies  those  integers  such 

1C  ^ 

that  Pk's  exceed  some  preselected  value  Pk.  We  let  N(h,v,y)  be  the 
number  of  integers  in  Ip^. 

A  /\ 

A  reasonable  defense  policy  is  to  set  M,  Pk  and  r.  Now 


M*(h,v,y) 

P£(h,vjy;r)  =  Y,  (1-(1-Pk(h,v,y,im))r)/M*(h,v,y) 
m=l 


where 


M*(h,v,y)  =  min(M(h,v,y) ,N(h,v,y)) 


so 


K*(h,v;r) 

P 


1 

R(h) 


R(h) 

M*(h,v,y)P*(h,v,y  ;r)dy 


0 


We  let 


M*(h, v) 


1 

R(h) 


R(h) 

M*(h,v,y)dy 


0 
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-* , 


and  define  P^(h,v;r), 


so 


P*(h,v,r)  =  M*(h,v)P*(h,v;r) 


P,(h,v;r) 
W*(h,v;r)  - - - 


where 


W*(h,v;r)  =  rM*(h,v) 


4.1.3  The  Effect  of  Penetration  Force  Size 

The  above  calculations  treat  the  penetration  force  as  a  point  which 
moves  through  the  region  of  defense  effectiveness.  For  a  small  number 
of  tightly  grouped  penetrators  this  may  be  adequate,  but  for  larger 
penetration  forces  practicing  jinking  and  maneuvering  and  trying  to 
control  decoys  the  size  of  the  force  may  be  an  important  consideration. 
Generally  speaking,  a  larger  penetrating  force  will  increase  the  firing 
opportunities  available  to  the  defense  and  perhaps  allow  him  to  select 
opportunities  with  higher  P^'s. 

The  easiest  way  to  account  for  this  effect  without  complicating  the 
model  is  to  assume  the  penetration  force  of  size  N  occupies  a  certain 
volume  in  space  which  does  not  change  as  penetrators  are  killed.  Thus, 
this  fixed  volume  implies  that  M(h,v,y)  and  Pk(h,v,y,m)  may  be  determined 
as  functions  of  N,  the  size  of  the  penetrating  force.  Regrouping  to 
present  a  smaller  volume  as  penetrators  are  killed  would  make  M  and  P^ 
depend  on  the  number  alive  at  any  instant  of  time.  This  complicates 
the  modeling  considerably  and  will  not  be  considered  in  this  report. 

4.1.4  Interceptor  Aircraft 

If  the  defense  site  is  an  interceptor  aircraft,  then  kill  potential 
may  be  computed  in  much  the  same  way.  We  set 

Kp(h,v)  =  M(h,v)Pk(h,v) 

for  aircraft  where  M(h,v)  is  the  number  of  passes  and  P^ChjV)  is  the 
kill  probability  of  a  pass.  We  allow  both  of  these  quantities  to 
depend  on  altitude  and  speed  of  the  penetration  force.  As  before,  we  can 
also  make  these  quantities  dependent  on  N,  which  may  increase  M  and 
perhaps  P^. 
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4.2  CUMULATIVE  KILL  POTENTIAL 

Since  kill  potential  over  a  number  of  defense  sites  is  additive, 
it  only  remains  to  determine  the  number  of  sites  encountered  in  a 
penetration  of  depth  x.  We  let  f(y)  be  the  density  of  defense  sites  at 
penetration  position  y,  then  the  number  of  sites  encountered  is 


N(x)  =  2R(h)  /  f(y)dy 


and  the  cumulative  kill  potential  is 


H(x)  =  2R(h)Kp(h,v)  J  f(y)dy 

0 


If  we  allow  h  and  v  to  depend  upon  penetrator  force  location. 


H(x)  =2  J  R(h(y))Kp(h(y) ,v(y))f(y)dy 


which  implies 


X(x)  =  2R(h(x))Kp(h(x) ,v(x))f(x) 


4.2.1  Nonhomogeneous  Defense 

Cumulative  kill  potential  is  easily  extended  for  a  mixture  of 
different  defense  types.  Once  again,  since  kill  potential  adds,  we 
merely  index  the  above  results  for  defense  type  k.  Thus 


H,  (x)  =  2 


A 


(h(y))K  (h(y)  ,v(y))f,(y)dy 
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with 


Ak(x)  =  2Rk(h(x))Kpk(h(x) ,v(x))fk(x) 


The  total  cumulative  kill  potential  is 
K 

H(x)  =  Kk(x) 
k=l 

which  implies  that 

K 

^(x)  =  ^ 

k=l 


4.3  ALLOCATION  OF  DEFENSE  CAPABILITY  TO  PENETRATORS 
OF  DIFFERENT  TYPES 

Penetrators  of  different  types  may  have  markedly  different  penetra 
tion  capabilities.  This  suggests  making  the  above  results  dependent  on 
penetrator  type.  Thus,  for  penetrator  type  j,  we  have 

Akj(x)  =  2^  (h(x))Kpkj(h(x)  ,v(x))fk(x) 

as  the  kill  intensity  of  type  k  defenses  against  type  j  penetrators. 
With 

-X 

V°  ■  J  Vy)dy 

0 

as  the  cumulative  kill  potential  if  type  k  defenses  are  used  against 
type  j  penetrators.  The  total  cumulative  kill  potential  that  may  be 
brought  to  bear  against  type  j  penetrators  is 
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H.j(l°  '£  Hkj(x) 


The  reader  should  compare  H(x)  in  4.2  and  H.j(x)  and  realize  that  H(x) 
is  H,j(x)  when  one  type  of  penetrator  is  considered. 

4.3.1  Fixed  Allocation  of  Defense  for  a  Nonhomogeneous 
Penetration  Force 

We  want  to  allow  the  penetration  force  to  be  made  up  of  penetrators 
of  different  types.  For  the  moment  let  us  assume  that  the  defense  fixes 
its  allocation  of  defense  to  penetrators  as  a  function  of  penetrator  force 
location.  We  let  (%j (x)  be  the  fraction  of  defense  capability  of  type  k 
allocated  to  penetrators  of  type  j  at  penetration  location  x,  where 


0  <  akj(x)  <  1 


all  k  and  j 


Z  aki(x)  s  1 


all  k 


Hkja(x)  ‘  /  Xkj<y)V|(y)dy 


H.ja(x)  ’£  /  Xkj(y)cl kj(y)dy 


k  0 

is  the  cumulative  kill  potential  allocated  to  penetrators  of  type  j 
The  kill  intensity  is 


I  I 


X-ja(x)  '  £  Xkj(x)tVj(x) 
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r 

I 


We  note  that  these  results  still  allow  one  to  use  the  models  in 
Section  3  which  require,  at  most,  the  numerical  integration  of 
A.ja(x)  or  oj.ja  and  4>.ja  functions  related  to  it.  That  is,  numerical 
integration  of  the  basic  differential  equations  and  problems  of  de¬ 
pendence  among  penetrator  types  have  still  been  avoided. 

4.3.2  Classification  of  Penetrator  Type  and  Adaptive 
Allocation  of  Defense 

If  one  can  distinguish  penetrators  of  different  types,  then  a^j(x) 
may  be  specified.  Classification  may  not  be  perfect;  moreover, 
a^j (x)  should  be  0  if  no  penetrators  of  type  j  are  alive. 


We  now  consider  an  allocation  process  which  takes  these  factors 
into  account.  We  assume  a  penetrator  of  type  j  is  classified  as  type  r 
with  probability  Pjr(x).  If  a  penetrator  of  type  r  has  relative  weight 
Wr(x)  and  if  ij  penetrators  of  type  j  are  alive  at  x,  then 


\j<x)  *  Vj  E  Pjr(x)Wr(x) 


where 


W*(x)  =  ^  Pjr(x)Wr(x) 


and  is  chosen,  so 


E  ■  ek(x)  5 1 


where  S^Cx)  is  a  defense  decision  function  relating  to  the  use  of  k  type 
weapons .  Thus , 


NWC  TP  5841 


i,W*(x)6,  (x) 

a  (x)  - — 

J  £  i4W*(x) 


j 


j  j 


°kj(x)  “  Ck^x;il,i2,'”’iJ)ijWj 


(x) 


where 


Ck(x;il»’”*1j)  = 


E 


-1 


L  J 


ek(x) 


In  this  section  we  find  ourselves  in  the  dependent  case  with 


K 

Aj  (x;i1,i2,  •  •  •  ,ij)  =  ijW*(x)  ^  A^  (xK^xjip  •  •  •  ,ij) 

k-1 


Also,  one  is  immediately  involved  in  the  numerical  solution  of  stochastic- 
differential  equations. 


4.3.3  Adaptive  Allocation  of  Defense  With 
Independence  of  Penetrator  Type 

We  have  seen  that  adaptive  allocation  of  defense  leads  to  dependence 
among  types  of  penetrator  groups.  Moreover,  to  carry  out  the  solution 
would  require  the  numerical  solution  of  the  appropriate  stochastic- 
differential  equations.  The  difficulty  generated  in  4.3.2  is  similar 
to  complexities  that  arise  in  the  analysis  of  penetration  aid  effects. 

We  nov,  consider  several  methods  of  handling  this  problem. 


I 


have 


The  fir, it  step  comes  in  analyzing  where  the  difficulty  occurs.  We 


K 


Aj(x;i1,***,xJ)  -  2  Vj(x)*kj(x) 


k=l 


I 
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where 


(x) 


2Rfc1(h(x))K  kj (h(x) ,v(x))fk(x) 
ijW*(x)6k(x) 

E  i.W^x) 

j 


The  problem  is  not  in  the  kill  intensity  portion  of  the  calculation, 
but  in  the  allocation  portion.  The  problem  is  really  one  of  approximating 
akj  (x)  . 

If  we  can  assume  that,  for  all  j,  there  is  an  attrition  function 
a(x) ,  independent  of  j  so  that 

ij  (x)  =  Nj  a(x) 

then 

N.W*(x)0,  (x) 

(x)  =  - - - 

E  Njw*(x) 

Now  Aj (x;ii , • •  •  ,ij)  =  Aj(x),  so  independence  holds.  In  fact,  the  model 
of  3.1  or  3.2  may  be  applied  to  yield  Qij (x)  for  each  j  independently. 

One  might  be  tempted  to  set 


‘Vj 


(x) 


i  W*(x)6k(x) 
E  NjW*(x) 


and  then  use  the  model  of  3.2.  This  is  certainly  possible,  but  the 
intent  of  model  3.2  is  to  consider  those  cases  where  A-j  (x)  may  be 
reduced  at  the  discretion  of  the  defense  when  penetration  force  size  is 
small.  Another  problem  is  that  the  above  expression  summed  on  j  does 
not  give  9k(x) . 


3t 


Another  possibility  is  to  replace  akj (x)  by  its  expected  value 


ECo^  (x)) 


\j 


(x)Q. 


1* 


Ll* 


Even  if  Q-j^  _  _  _  ^j(x)  were  known,  this  procedure  would  suffer  from  the 
fact  that,  in  general, 


E(\j  (x))  ^  0ku) 


If,  however,  we  define  Ej (x)  as  the  expected  number  of  survivors 
of  type  j  alive  at  x  and  set 


E  (x)W*(x)6k(x) 

E  E  (x)W*(x) 
j 


we  have 


Z  ■  ekw 

j 

The  only  problem  remaining  is  that  Ej  (x)  is  not  known. 

We  turn  to  this  problem  under  more  general  conditions  in  5.1,  but 
for  now  consider  setting 


N.W*(x)0.  (x) 

J  3  k 

N.W*(x) 

2  2 


Using  the  model  of  3.1,  3.2,  or  3.4,  whichever  is  appropriate,  we  obtain 
(x)  and  compute 

^  J 
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Ej1)(x) 


-  E 


/n  (x)W*(x)0,  (x) 

5*7  00  «  vvv3 — - — 

j  £Ej(1)(x)Wj(x) 

(2)  d') 

which  allows  one  to  obtain  Q.  (x)  and  hence  (x) .  This  procedure  is 

repeated  until  changes  in  the  Ej (x)  calculation  are  acceptable.  This 
procedure  will  be  expanded  upon  in  Section  5.1. 


Section  5.  PENETRATION  AIDS  AND  TACTICS 


We  now  consider  the  effect  of  various  penetration  aids  and  tactics 
on  the  kill  potential  the  defense  can  mount.  We  shall  consider  trajectory 
variations,  decoys,  ECM,  defense  suppression,  self-defense,  and  standoff 
tactics. 

The  reader  should  recall  that,  with  reasonable  assumptions,  we  have 
been  able  to  maintain  independence  of  penetration  types.  As  pointed 
out  in  Section  3,  if  one  is  forced  to  go  to  dependence,  then  the  computing 
may  become  excessive  for  large  penetrations.  In  this  section  it  is 
virtually  impossible  to  avoid  dependence  because  of  the  synergistic 
effects  of  one  type  of  penetration  group  with  another.  Indeed  these 
effects  are  an  important  part  of  penetration  aid  analysis.  Before 
treating  the  penetration  aid  problem  we  consider  a  general  approximation 
method  which  maintains  independence  while  accounting  for  synergistic 
effects  in  an  average  value  way. 

5.1  RECURSIVE  SOLUTION  PROCEDURE  TO  MAINTAIN  INDEPENDENCE 
AND  ALLOW  ASSESSMENT  OF  SYNERGISTIC  EFFECTS 

Penetration  aid  analvsis  will  lead  to  Aj  which  depend  not  just  on 
ij  but  also  on  ii>i2»'**»  anc*  ij.  We  recall,  however,  that  cumulative 
expected  kill  potential  has  always  been  the  parameter  of  interest  to  us. 
Thus,  taking  Aj  (x;ij_,  •  •  •  ,ij)  and  replacing  it  by 


E[Aj (x) ] 


presents  no  conceptual  difficulty.  The  problem  is  that  Qi^ , . . . ,ij(x) 
is  unknown. 

Assuming  that  one  has  an  initial  estimate  E[A^(x)]  of  E[A.«(x)] 

J  (0)  J 

then  by  using  a  model  of  3.1,  3.2,  or  3.4  one  may  obtain  (x) ,  since 
independence  holds.  In  general,  at  the  kth  iteration  one  obtains 


E[A<k)(x)]  =  V  A.(x;i1,...,iT)Qifk)(x),---,Q.(k)(x) 
J  J  1  J  11  \j 


1  •  •  •  i 


.(k+1) 


from  which  Q.  (x)  is  obtained  using  a  model  of  3.1,  3.2,  or  3.4. 


This  iterative  procedure  is  continued  until 


max 

x,j 


(k) 


X<k+1> (x) 


<  e 


where  e  is  a  small  positive  number. 

Another  possibility  is  to  approximate  Aj  by  Aj  (x;Ej_(x)  ,  •  •  •  ,Ej(x)) 
where  Ej (x)  is  the  expected  number  of  survivors  of  type  j  alive  at  x. 

In  this  case  we  begin  with  an  approximate  solution  Ej(^)(x)  which  yields 
Qij (D (x) .  At  the  kth  iteration,  Aj (x) , • • • (x) )  yields 
Qij(k+l)(x),  using  a  model  of  3.1,  3.2,  or  3. A.  Now 

E<t+1)(x)  -Ji.qWw 

J  J  j 


This  iterative  procedure  stops  when 


max 

j 


E(k) 

J 


(x)  -  E 


(k+1) 

j 


<  £ 


In  this  report  we  shall  concentrate  on  this  latter  approach.  We 
have  already  seen  in  4.3.3  that  it  produces  less  complexity  than  dealing 
with  E(Aj(x)).  It  still  satisfies  the  requirement  of  allowing  for  the 
evaluation  of  synergistic  effects  among  penetration  aids. 


5.2  TRAJECTORY  VARIATION 

The  model  already  allows  altitude  h  and  speed  v  to  depend  on 
penetrator  location  x.  If  one  has  available  a  model  for  the  effects 
of  jinking  and/or  terrain  following,  inputs  to  this  model  in  terms  of 
M  and  for  a  single  defense  site  may  be  modified  accordingly.  If 
M  and/or  are  not  reduced  sufficiently,  then  the  effect  of  maneuvering 
may  not  reduce  the  kill  potential  of  the  site  below  the  level  the  defense 
intends  to  expend  in  the  first  place.  In  addition,  3,  the  reliability 
failure  rate,  will  increase  under  a  maneuvering  tactic  and  the  average 
speed  of  penetration  will  decrease,  further  increasing  y  =  3/v. 

Thus  the  model  as  structured  allows  for  jinking  and  terrain¬ 
following  inputs  from  other  sources  and  allows  for  an  evaluation  of 
these  tactics  through  their  effects  on  M,  Pfc,  and  y,  the  reliability 
factor. 
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5.3  DECOYS 

The  model  easily  handles  decoys  by  allowing  one  of  the  penetrating 
subgroups  to  be  a  group  of  decoys.  Important  parameters  in  the  model 
with  regard  to  the  decoys  subgroup  j  are  Uj ,  the  reliability  factor; 
Wj(x),  the  weight  the  defense  places  on  decoys;  and  Pjk(x) ,  the  proba¬ 
bility  that  a  decoy  is  classified  as  a  penetrator  of  type  k.  It  is  not 
at  all  obvious  that  the  defense  should  place  a  weight  of  0  on  decoys. 

In  fact,  the  evaluation  of  defense  policy  with  regard  to  assignment 
of  all  weights  is  an  important  attribute  of  the  model. 

5.4  ELECTRONIC  COUNTERMEASURES 

Electronic  countermeasures  (ECM)  may  affect  defense  kill  potential 
in  several  ways.  The  user  may  specify  any  or  all  of  the  following 
possibilities . 

5.4.1  Noise  Jamming 


If  there  are  ij  penetrators  of  type  j  alive  at  x,  P(j)  is  the  jamming 
power  of  a  penetrator  of  type  j,  and  A(j)  is  the  average  radar  reflective 
area  of  a  type  j  penetrator,  then  we  define  expected  noise  jamming  power 
of  the  group  at  x  to  be 

J 

P(x)  =  £  Ej(x)P(j) 

j=l 


With  simplifying  assumptions  regarding  penetrator  volume  and  resolvability 
of  the  defense  radar,  the  expected  radar  reflective  area  of  the  group  at  x  is 

J 

A(x)  =  £  Ej(x)A(j) 

J-l 

where  Ej(x)  is  the  expected  number  of  penetrators  of  type  j  alive  at  x. 

Essentially  noise  jamming  reduces  the  range  at  which  the  defense 
can  take  action.  If  Rk(x)  is  the  burnthrough  range  at  x  for  defense 
type  k,  then  we  estimate2  R^(x)  by 


2  Naval  Operations  Analysis.  U.  S.  Naval  Institute,  Annapolis,  Md. , 
1968.  (Library  of  Congress  Card  67-28645.) 
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VE< 


x))  =  / 


PT(k)G(k)A(x) 


4ttP(x) 


where  PfOO  is  the  transmitted  power  associated  with  a  defense  of  type,  k, 
G(k)  is  the  associated  antenna  gain,  ^nd 


E(x)  =  [E1(x) , • • ' >Ej(x) ] 


Thus  the  effective  range  of  a  defense  of  type  k  against  a  penetrator 
of  type  j  is 


Rk..(h(x),E)  =  maxjO,min[Rk(E(x))  -  vAT^  .R^  (h(x)  )  ] 


where  ATj^j  is  the  delay  for  the  system  to  acquire  once  burnthjough  has 
occurred.  Now  in  computing  kill  intensity  one  uses  R]cj(h(x),E) 
in  place  of  Rkj(h(x)).  If  one  can  reasonably  assume  that  jamming 
reduces  the  basic  single-shot  Pj^s,  this  effect  can  also  be  considered. 


Thus  with  noise  jamming 


Akj(x,E)  =  2Rkj(h(x),E)Kpkj(h(x),v(x),E)fk(x) 


where 


\  .(h(x),E) 

r  K.J 


K  ,  (h(x) ,v(x) ,E)  = 


Pkj 


M.  .  (h(x)  , v(x)  ,E,y)P^^  (h(x)  ,v(x)  ,E,y)dy 
■W.W  n  3 


(h(«/,ij/  q 


Kk  j 


and  the  expressions  under  the  integral  sign  are  obtained  using  Rkj(h(x),E) 
We  also  have  Kpkj  (h(x)  ,  v(x)  ,?)  =  M^j  (h(x)  ,v(x)  ,E)PKkj  (h(x)  ,  v(x)  ,E)  where 


R,  (h(x)  ,£) 


"kJ 


(h(x) ,v(x) ,E)  = 


\j (h(x>  »E>  0 


(h(x) ,v(x) ,E,y)dy 
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and  P^j  (h(x)  ,v(x)  ,E)  is  determined  accordingly. 

Even  though  there  is  an  obvious  cancellation  of  R^j  terms ,  it  is 
instructive  to  maintain  the  form 


Akj(x,E)  -  2Rkj(h(x)  .E)^  (h(x)  ,v(x)  ,l)PKkj  (h(x)  ,v(x)  ,E)fk(x) 

and  note  that  noise  jamming  essentially  degrades  the  range  of  effectiveness 
of  a  site,  thus  changing  M  and  P^  for  a  single  site  and  reducing  the 
number  of  sites  encountered. 

5.4.2  Deception 


One  important  aspect  of  ECM  is  to  make  the  classification  problem 
of  the  defense  more  difficult.  If  a  particular  class,  j*,  is  devoted 
to  providing  an  ECM  capability,  then  we  may  allow  the  classification 
probability  Pjr(x)  to  depend  on  ij*,  the  number  in  this  subgroup. 
Referring  to  Section  4.3.3,  we  have 

iiWj(x;ij*)6k(x) 


Okj(x)  - 


where 


w*(x;ij*)  =  S  Pjr(x;ij*)wr(x) 


Proceeding  as  in  4.3.3,  we  have 


where 


(x;EjA)  = 


E  (x)W*(x;E  (x))6,(x) 

J  3 _ -  J_ _ £ - 

J 

(x)W:  ;E  .*(x)) 

j=l  J  J 


£  PjrvA;Ej*(x))Wr(x) 
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Now  we  define 


yielding  kill  intensity  due  to  jamming,  deception,  and  defensive 
allocation. 

There  will  be  an  attempt  throughout  the  remainder  of  the  report  to 
make  the  notation  as  descriptive  as  possible  of  the  effects  that  are 
being  considered.  Thus 

xkJa(x3)  -  5y(x)Xkj(x>S) 

describes  a  model  ^n  which  the  kill  intensity  depends  on  allocation  (a) 
and  noise  jamming  E,  but  deception  due  to  Ej*  is  not  considered. 

5.5  DEFENSE  SUPPRESSION 

We  now  want  to  consider  the  possibility  of  using  weapons  to  suppress 
or  kill  defense  sites.  We  assume  each  penetrator  of  type  j  carries 
Sjk  weapons  to  either  kill  or  suppress  defense  units  of  type  k.  With 
an  appropriate  definition  of  penetrator  subgroup  of  type  j  ,  it  is  assumed 
that  each  (j  ,k)  combination  may  be  categorized  as  "suppression"  or 
"destruction."  Thus,  a  given  penetrator  type  j  will  not  both  suppress 
and  destroy  a  defense  of  given  type. 

We  assume  that  these  weapons  are  expended  at  a  rate  proportional  to 
fk(x) ,  the  density  of  type  k  defenses  at  penetratxon  depth  x.  Moreover, 
it  is  assumed  that  these  weapons  are  expended  over  (0,Zk).  Should  one 
want  to  use  these  weapons  over  (Yk,Zk)  where  0  <  Yk  <  Zk  then  we  replace 
k  type  defenses  by  two  identical  defenses  k^  and  V-2  where 


fk  (x)  =  fk(x) 

o  <  X  <  Y] 

l 

=  0 

Y.  <  x 
k 

fk  (x)  =  0 

0  <  x  <  Y. 

2 

/ — > 
x 

II 

Y.  <  x 
k 

'sssasammm 
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Now  Sjj^  =  0 
Since  f^Cx) 
allocated  to 


and  SjjC2  =  sjk  with  the  weapons  allocated  over  (0,Zk)  . 
=  0  for  0  <  x  <  Yk,  the  Sjk  weapons  per  aircraft  are 
k  type  defenses  over  (Yk,Zk)  . 


Clearly  a  suppression  weapon  must  be  used  at  some  point  before  x 
to  be  effective  at  x.  We  shall  make  the  simplifying  assumptions  that  the 
number  of  penetrators  alive  at  x  approximates  those  alive  at  a  position 
to  launch  suppression  weapons  at  defenses  at  x.  Based  on  the  expected 
number  of  surviving  defense  sites,  fk(x,S),  the  density  of  sites  after 
suppression,  satisfies 


2 


/ 


Rk(h(y))fk(y,S)dy 


)fk(y)dy 


r(y)Nk<zk) 


Ej(y)Rk(h(y))fk(y)dy 


where  r(y)  is  the  smallest  integer  such  that  the  expression  on  the  right 
is  non-negative,  Rk(h)  is  a  nominal  radius  of  action  for  the  attackers 
against  defense  type  k  with  the  constraint 

\(h(y))  *  m^X  ^.(hCy)) 


So  that  the  attackers  attempt  to  suppress  all  sites  that  can  have  an 
effect  on  the  penetration 

A 

Nk(Zk)  =  2  J  Rk(h(y))fk(y)dy 
0 


This  result  follows  from  the  assumption  that  the  rate  of  weapon  expenditure 
is  proportional  to  fk(x)  and  the  requirement  that  all  weapons  be  expended 
by  the  time  Zk  is  reached.  The  smaller  Rk(h(y)),  the  more  effective  the 
penetration.  Thus,  if  ECM  is  utilized 

Rk(h(x))  >  mJX  Rk.(x,l) 
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and  defense  suppression  together  with  ECM  can  have  a  beneficial  joint 
effect. 

Taking  the  derivative  with  respect  to  x  yields 


fk(x,S)  = 


fk(x) 


-£ 


(1‘<1‘psjk> 


r(x) 


)ej<*> 


r(x)Nk(Zk) 


0  <  x  <  Z, 


=  fk(x) 


Zk  S  x 


where  r(x)  is  the  smallest  integer  such  that  the  expression  in  brackets 
is  positive.  In  most  instances  we  can  set  r  =  1  and  then 


fk(x,S)  =  max 


0,fk(x) 


i .  y  .jZsikgi [X1 

~  \(zk> 


0  <  x  £  Z, 


=  fk(x) 


Zk^x 


Now 


Akj(x,E,S)  =  2Rkj(h(x),E)Mkj(h(x),v(x),l) 
X  ^Kkj(h(x)»v(x) >E)fk(x,S) 


and 


^kjot^x*^*^  -  *^kj  (x ’ *^kj  ^x 
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5.6  SELF-DEFENSE 

The  situation  in  self-defense  is  much  the  same  as  in  defense  sup¬ 
pression.  In  this  case  the  weapons  that  j  type  penetrators  allocate  for 
defense  against  the  weapons  that  k  type  defenses  utilize  must  be  used 
against  these  types  of  weapons. 

If  Wj,j(h,v)  is  the  expected  number  of  weapons  a  type  k  defense  can 
get  off  at  a  type  j  penetrator,  then  the  total  number  of  weapons  expended 
by  k  type  defenses  against  j  type  penetrators  is 

Wkja('X’^,Ej*’S^  \j(h(y)  ,v(y)  ,1)11^  (h(y)  ,E)a^  (y,E^)fk(y  ,S)dy 

0 


if  jamming,  deception,  suppression  and  defense  allocation  are  all 
considered. 

The  penetrators  of  type  j  attempt  to  counter  these  weapons  with 
djk  defense  weapons  per  j  type  penetrator  with  kill  probability  P^jk. 
per  defense  weapon.  We  assume  that  the  j  type  penetrators  expend  these 
weapons  at  a  rate  proportional  to  the  density  of  weapons  defense  type  k 
launches  at  penetrator  type  j.  Thus,  these  weapons  are  expended  at  a 
rate  proportional  to  dWkjo[(x,E,Ej*,S) /dx  over  distance  (0,Zjk)  .  The 
cumulative  number  of  surviving  weapons  due  to  self-defense  is 


Wkja^x’E,Ej*,S,D^  ~  2  /  (b(y)  ,v(y)  ,E  ,D)  (h(y)  ,1)^  (y  ,Ejsl()  fR(y  ,S)  dy 

0 


/ 


max 


0, 


/  _  d1kPdikEl(y) 


dy 


w 


j 


dy 


m  anurmin  k  ton*  StmaaUmt&illm 
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It  follows  that 


dW  (x,E,E  S,D)  r  dfc.  .  (x, i,E44,S> 


■^.....j.......  rkja(x  E,Ej<,S) 

dx  L  dx  \  VV 


where 


W  ■ 


=  W,  .  (Z.,  ,E,E.  .  ,S) 
kjor  jk’  *  2* 


It  follows  that 


->  *  /  d  .  PE. 

M^OKxJ.vCxM.D)  "maxO.^^.vCxJ.E)  ( 1  ~  N  (Z  J 

\  jk'  jk' 


0  <  x  <  Z 


Mkj(h(x),v(x),E) 


Z <  x 
jk 


Akj(x,E,S,D)  =  2Rkj(h(x),S)Mkj(h(x),v(x),E,D) 


X  PKkj fV^x^ »E)fk(*>S) 


Xkja(x’P’S’D)  =  \j(x,E,S,D)akj(x,E,E  ) 
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5.7  STANDOFF  WEAPONS 

We  assume  that  one  particular  subgroup,  jg,  say,  intends  to  launch 
standoff  weapons  at  Dg  for  impact  at  rg  >  Dg.  If  each  element  of 
subgroup  jg  carries  wjq  standoff  weapons,  then  the  probability  that 

Wj0ij0  standoff  weapons  are  launched  at  Dg  is  Qij^(Dg).  The  probability 

distribution  of  the  number  of  standoff  weapons  surviving  to  rg  is 


j0 

Vro>  "  £ 

ij0“[i/uJ, 


Qi.  (D0)Qi|w.  i  (D0,r0) 
*  Jg  Jg  Jg 


where  |i/wjQJ*  is  the  smallest  integer  >  i/wjQ. 
3.1  may  be  used, 


1|W1  1i 
J0  J0 


(Dg » rg) 


[H(D0,r0)] 


(w  i  -  i)  ! 
J0  J0 


Providing  the  model  of 


' 

-H(D0,rg) 

e 


wJoxjo 

1  '  2  Qi|w.  i.  (D0,r0) 

i=l  Jo  Jo 


H<Do*ro) 


Vy)dy 


where  Am(y)  is  the  kill  intensity  due  to  all  defenses  against  standoff 
missiles  at  position  y.  It  is  understood  that  Am(y)  is  obtained  using 
concepts  developed  in  Section  4.0. 
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Appendix  A 

A  NUMERICAL  COMPARISON  OF  POISSON  AND 
BINOMIAL  DISTRIBUTIONS 


In  this  appendix  we  compare  the  truncated  Poisson  with  the  binomial 
distributions  associated  with  SLS  and  MSE  firing  doctrines.  For  a  fixed 
H,  there  are  a  number  of  combinations  of  M  and  P^  that  may  be  used.  The 
binomial  results  depend  on  both  M  and  Pfc  and  it  is  important  to  see  how 
the  truncated  Poisson  results  compare  with  more  acceptable  binomial  results 

Tables  A-l  through  A-7  allow  one  to  make  this  comparison  for  a 
range  of  N,M,Pk  and  H.  Generally  speaking  the  truncated  Poisson  results 
lie  between  those  for  SLS  and  those  for  MSE.  Since  most  penetrations 
will  involve  a  mixture  of  MSE  and  SLS  firing  doctrines,  this  is  an 
encouraging  result.  Stated  differently,  firing  doctrine  can  produce 
a  larger  difference  in  results  for  binomial  models  than  the  difference 
between  a  truncated  Poisson  and  either  binomial  model. 


TABLE  A-l.  Comparison  of  Binomials  and  Truncated  Poisson  for  H  =  0.1 


Expected  Kills  for  the  Binomials  SLS/MSE 


.005  (20) 

■  095/ ( -095) 

. 1 00/ ( .098) 

. 100/ ( .099) 

.ioo/Moo) 

.01  (10) 

.096/(.096) 

. 100/( .098) 

. )00/( . 100) 

,I00/(.100) 

■025  W 

•096/(.096) 

. ico/( .099) 

.  ioo/( . 100) 

■I00/(.I00) 

•05  (2) 

. 098/  ( .098) 

. 100/( . 100) 

.  ioo/Moo) 

.ioo/Moo) 

•10  (1) 

. 100/( . 100) 

. !00/( . 100) 

_ 

.  100/( . 100) 

.I00/(.I00) 

Properties  of  Truncated  Poisson 


Expected 
number 
killed 

For  SLS: 

Max  $  error 
Avg  X  error 

For  MSE: 

Max  $  error 
Avg  $  error 

Avg  $  error 
for  both 
SLS  and  MSE 


IOO/( . 100) 

ioo/Moo) 
100/( . 100) 
IOO/( . 100) 
!00/( . 100) 
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TABLE  A-2.  Comparison  of  Binomial',  and  Truncated  Poisson  for  H  =  0.25. 


Expected  Kills  for  the  Binomials  SLS/MSE 


.248/ (.236 
,248/( .236 


. 25/ ( .244) 
•25/( .245) 
.25/ ( .247) 
■25/( .25) 
.25/( .25) 

.  25/ ( . 25) 


.  25/ ( .248) 
,25/( .248) 
25/(.25) 

,  25/( .25) 
25/( .25) 

.  25/ ( . 25) 


25/( 

.249) 

25/ ( 

.25) 

25/( 

.25) 

25/ ( 

.25) 

25/ ( 

25) 

25/(. 

.25) 

Properties  of  Truncated  Poisson 


Expected 
number 
ki  lied 


For  SLS: 

Max  I  error 
Avg  $  error 

For  MSE: 

Max  $  error 
Avg  $  error 

Avg  $  error 
for  both 
SLS  and  MSE 


TABLE  A-3.  Comparison  of  Binomials  and  Truncated  Poisson  for  H  =  0.50. 


.005  (100) 
.01  (50) 
.025  (20) 
.05  (10) 

.10  (5) 

■25  (2) 

•50  (I) 


Expected  Kills  for  the  Binomials  SLS/MSE 


.5/(.48g) 
.5/ (.490) 
•5/( .494) 
-  5/ ( • 5) 
■5/ (-5) 

. 5/ ( • 5) 
•5/( -5) 


•394/(.394) 

.484/( .443) 

.  5/ ( .477) 

. 395/ ( -  395) 

. 484/ ( .444) 

•5/(.478) 

.3)7/ (-397) 

.486/( .447) 

•5/(.482) 

. 401/ ( .401) 

.487/ (.450) 

•5/( .49) 

. 4 1 / ( .41 ) 

•49/(.46) 

•  5/(  .5) 

•437/ ( .437) 

. 50/ ( . 50) 

•5/ (.5) 

•5/ (.5) 

•5/( -5) 

•  5/ (.5) 

. 5/( • 495) 
.5/ (.496) 
. 5/ ( -  5) 

-  5/ ( -  5) 

•  5/ ( -  5) 
•5/ (.5) 
•5/( .5) 


Properties  of  Truncated  Poisson 


Expected 
number 
ki  lied 


For  SLS: 

Max  %  error 
Avg  %  error 

For  MSE: 

Max  %  error 
Avg  %  error 

Avg  I  error 
for  both 
SLS  and  MSE 


•  393 

.484 

.5 

21 .4 

3.2 

0 

5.6 

1  .2 

0 

21 .4 

-9.2 

-4.8 

5.6 

-4.6 

-2.2 

5.6 

-1.7 

-1.1 

mm 


mvif: M  m 
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TABLE  A4.  Comparison  of  Binomials  and  Truncated  Poisson  for  H  =  1. 


N 

pk(mT^^ 

1 

! - 

2 

5 

10 

20 

Expected 

Kills  for  the 

Binomials  SLS/MSE 

.005  (200) 

•633/ ( .633) 

.897/ ( . 788) 

•999/( .908) 

1  .00/ ( . 954 ) 

1  .oo/( .978) 

.01  (100) 

.634/( .634) 

.898/ (.790) 

•999/ ( .910) 

1 . 00/ ( .956) 

1  .oo/( .98) 

.025  (40) 

.637/( -637) 

. 90 1 /( -795) 

1 . 0/ (.917) 

1  .o/( .963) 

1 .o/( .987) 

.05  (20) 

.642/( .642) 

•906/( .802) 

1 .o/( .927) 

1  .o/(  .975) 

1  .o/O.o) 

.10  (10) 

.65 1 /( .651) 

■  9 1 5/ (.819) 

1 ■ 0/ ( .95) 

1. 0/(1 .0) 

1. 0/0.0) 

•  25  (<•) 

. 684/ ( .684) 

. 945/ ( -  875) 

l .0/(1 .0) 

i.o/(i.o) 

1. o/O.o) 

•50  (2) 

•  75/ ( .75) 

1 .0/(1 .0) 

1 .0/(1 .0) 

1  .0/(1 .0) 

1. o/O.o) 

Properties  of  Truncated  Poisson 

Expected 

number 

killed 

.632 

.896 

.999 

1.0 

1  .0 

For  SIS: 

Max  %  error 

36.8 

I0.lt 

0.  1 

0 

0 

Avg  $  error 

8.2 

3.8 

0. 1 

0 

0 

For  MSE: 

Max  $  error 

36.8 

-13.7 

-10.0 

-4.8 

-2.2 

Avg  %  error 

8.2 

-5.4 

-5.1 

-2.0 

-0.2 

Avg  $  error 
for  both 

SLS  and  MSE 

8.2 

-0.8 

-2.5 

-1.0 

-0. 1 

TABLE  A-5.  Comparison  of  Binomials  and  Truncated  Poisson  for  H  =  2. 


N 

pK(Mr\^ 

1 

2 

5 

10 

20 

Expectec 

Kills  for  the 

Binomials  SLS/MSE 

•05  (40) 

•  87/087) 

1 .47/(1 .28) 

1.98/(1.68) 

2.00/(1.86) 

2. 00/(1 .95) 

.10  (20) 

.88/ (.88) 

1 .49/(1 .30) 

1.99  /  0.72) 

2. 00/(1 .90) 

2.00/(2.00) 

•25  (8) 

.90/090) 

1  .53/0.37) 

2. 00/(1 . 84) 

2.00/(2.00) 

2.00/(2.00) 

•  50  (<t) 

.94/< .94) 

1 .62/(1 .50) 

2.00/(2.00) 

2.00/(2.00) 

2.00/(2.00) 

•6?  (3) 

■  96/ ( .96) 

1.70/0.62) 

2. 0/(2. 0) 

2.00/(2.00) 

2.00/(2.00) 

1.0  (2) 

1 .00/(1.00) 

2.00/(2.00) 

2. 0/(2. 0) 

2.00/(2.00) 

2.00/(2.00) 

Properties  of  Truncated  Poisson 

Expected 

number 

killed 

.87 

1 .46 

1  .98 

2.00 

2.00 

For  SLS: 

Max  %  error 

13.0 

27.0 

1 .0 

0 

0 

Avg  %  error 

KM 

9-7 

0.8 

0 

0 

For  MSE: 

Max  %  error 

27.O 

-17.8 

-7.5 

-2.6 

Avg  %  error 

1  ill * 

1.1 

-6.3 

-2.1 

-0.1( 

Avg  $  error 
for  both 

SLS  and  MSE 

H 

5. it 

-2.7 

-1.1 

-0.2 

1 
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TABLE  A-6. 


\  N 


Comparison  of  Binomials  and  Truncated  Poisson  for  H  =  5. 


5 


.05  (100) 

•  99/ ( .99) 

•10  (50) 

1  .00/(1.00) 

.25  (20) 

1.00/(1.00) 

.50  (10) 

1 .00/(1.00) 

.67  (7.5) 

1 .00/(1.00) 

.83  (6) 

1.00/(1.00) 

1.00  (5) 

1.00/(1.00) 

Expected  Kills  for  the  Binomials  SLS/MSE 


^ . 98/(4 .01) 
^  .99/(4 . 10) 
4.99/(4.38) 
5.00/(5.00) 
5.00/(5.00) 
5.00/(5.00) 
5.00/(5.00) 


Properties  of  Truncated  Poisson 


Expected 

number 

kilted 

For  SLS: 

Max  $  error 
Avg  $  error 

For  MSE: 

Max  $  error 
Avg  $  error 

Avg  $  error 
for  both 
SLS  and  MSE 


5/(4. 52) 
5/(4 .63) 


5/(5. 0) 
5/(5. 0) 
5/(5. 0) 


.99 

1.95 

4.12 

4.98 

1.0 

2.5 

17.6 

0.4 

1.0 

1.8 

6.8 

0.3 

1  .0 

-5-4 

-28.3 

-24.2 

1  .0 

-1.2 

-9.0 

-8.2 

1  .0 

0.3 

-l.l 

-4.0 

TABLE  A-7.  Comparison  of  Binomials  and  Truncated  Poisson  for  H  =  10. 


10/(8.02) 

10/(8.19) 

10/(8.75) 


1 

2 

5 

10 

Expected 

Kills  for  the 

Binomials  SLS/MSE 

.05  (200) 

1 .00/(1 .00) 

2/(1. 99) 

4.96/(4.36) 

8.78/(6.42) 

.10  (100) 

1 .00/(1.00) 

2/(1. 99) 

4.97/(4.39) 

8.81/(6.51) 

•25  (40) 

1 .00/(1.00) 

2/0.99) 

4.98/(4.50) 

8.92/(6.84) 

.50  (20) 

1 .00/(1.00) 

2/0.99) 

4.99/(4.69) 

9.12/(7.50) 

•67  (15) 

1 .00/(1.00) 

2/(2. 00) 

5.00/(4.82) 

9.28/(8.08) 

■83  (12) 

1 .00/(1 .00) 

2/(2. 00) 

5.00/(4.93) 

9. 89/(8  84) 

1.00  (10) 

1 .00/(1 .00) 

2/(2. 00) 

5.00/(5.00) 

i0. 0/(10. 0) 

Properties  of  Truncated  Poissor 


Expected 

number 

killed 

For  SLS: 

Max  $  error 
Avg  %  error 

For  MSE: 

Max  3  error 
Avg  I  error 

Avg  I  error 
for  both 
SLS  and  MSE 
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Tables  A-8  through  A-14  were  constructed  by  considering  that  the 
results  of  any  penetration  are  indeed  a  mixture  of  MSE  and  SLS  results. 
Generally  speaking,  when  M/N  is  large  one  would  expect  SLS  to  predominate. 
We  consider  a  weighting  of  che  form 


4>(a)  =  e_aM/N  E(MSE)  +  (l  -  e"0*1^)  e(SLS) 

where  a  is  chosen  so  that  e  =  1/2  when  M/N  is  as  close  to  2  as  is 

possible.  Roughly,  we  assume  MSE  is  as  likely  as  SLS  when  M/N  =  2. 

Other  weightings  could  be  chosen,  but  the  point  is  to  compare  the 
results  of  the  truncated  Poisson  with  some  reasonable  combination  of 
SLS  and  MSE  results. 

The  truncated  Poisson  gives  results  fully  comparable  with  classical 
binominal  results  when  the  uncertainty  of  firing  doctrine  is  taken  into 
account.  The  approximation  is  worse  when 

1/2  N  <  H  <  2N 


TABLE  A-8.  Comparison  of  Binomials  and  Truncated  Poisson  for  H  =  0.1. 


N 

1 

2 

5 

10 

20 

Expected  Kills  for  the  Weighted  Binomials  SLS/MSE 


.005  (20) 

.095 

.100 

.100 

.  too 

.  100 

.01  (10) 

.096 

.100 

.  100 

.100 

.100 

.025  W 

.096 

.100 

.100 

.100 

.100 

•05  (2) 

.098 

.100 

.100 

.  100 

.  100 

•10  (1) 

.100 

.100 

.100 

.  too 

.100 

Properties  of  Truncated  Poisson 


Expected 

.095 

.100 

.100 

.100 

.100 

number 

ki lied 

Max  %  error 

5.0  at  0. 1 

0 

0 

0 

0 

Max  for  P^  <  0.5 

5.0  at  0. 1 

0 

0 

0 

0 

Min  error 

0  at  0.005 

0 

0 

0 

0 

Min  for  P«  <  0.5 

0  at  0.005 

0 

0 

0 

0 

Avg  %  error 

2.0 

0 

0 

0 

0 

PK  S  0.5 

* 

N 

pk(m)"-^^ 

1 

2 

5 

to 

20 

Expected  Kills  for  the  Weighted  Binomials 

SLS/HSE 

.005  (50) 

.222 

■11 

.250 

.250 

.25 

.01  (25) 

.222 

.249 

.249 

.25 

.025  (10) 

.221) 

.246 

.248 

.25 

.25 

•05  (5) 

.226 

.245 

•  25 

.25 

.25 

.10  (2.5) 

.232 

.248 

.25 

.25 

.25 

■25  (1) 

.250 

.250 

.25 

.25 

.25 

Properties  of  Truncated  Poisson 


Expected 

.221 

.248 

.25 

.25 

.25 

number 

killed 

Max  %  error 

11.6  at  0.25 

0.8  at  0.25 

0 

0 

0 

Max  for  Pk  <  0.5 

11  .6  at  0.25 

0.8  at  0.25 

0 

0 

0 

Min  %  error 

0.4  at  0.005 

-1.2  at  0.05 

-0.8  at  0.025 

-0.4  at  0.01 

0 

Min  for  Pk  <  0.5 

0.4  at  0.005 

-1.2  at  0.05 

-0.8  at  0.025 

-0.4  at  0.01 

0 

Avg  %  error 

3.5 

-0.2 

-0.2 

-0.  1 

0 

PK  4  0.5 

TABLE  A-10.  Comparison  of  Binomials  and  Truncated  Poisson  for  H  =  0.5. 


" . N 

pk(mP'\^ 

1 

2 

5 

10 

20 

Expected  Kills  for  the  Weighted  Binomials 

SLS/MSE 

.005  (100) 

.394 

.484 

.5 

.5 

.499 

.01  (50) 

.395 

.484 

.499 

.498 

.498 

.025  (20) 

.397 

.484 

.496 

.497 

.5 

.05  (10) 

.401 

.478 

.495 

.5 

.5 

.10  (5) 

.410 

.475 

.5 

.5 

.5 

■25  (2) 

.437 

.5 

.5 

.5 

.5 

•50  (1) 

.5 

.5 

.5 

.5 

.5 

Properties  of  Truncated  Poisson 


- 1 

Expected 

.393 

.484 

.5 

- , 

.5 

r- 

•  5 

number 

ki lied 

Max  %  error 

21.4  at  0.5 

3.2  at  0.5 

0 

0 

0 

Max  for  P|(  4  0.5 

21.4  at  0.5 

3.2  at  0.5 

0 

0 

0 

lii n  %  error 

0.2  at  0.005 

-1.9  at  0. 10 

-1.0  at  0.05 

-0.6  at  0.025 

-0.4  at  0.01 

Mi n  for  Pk  <  0.5 

0.2  at  0.005 

-1.9  at  0. 10 

-1.0  at  0.05 

-0.6  at  0.025 

-0.4  at  0.01 

Avg  %  error 

5-6 

0.5 

-0.3 

-0.1 

-0. 1 

Pk  4  5 
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TABLE  A-ll.  Comparison  of  Binomials  and  Truncated  Poisson  for  H  =  1. 


Expected  Kills  for  the  Weighted  Binomials  SLS/MSE 


.005  (200) 

.633 

.897 

.999 

1 .000 

.999 

.01  (100) 

.636 

.898 

.999 

•  999 

.996 

.025  (60) 

.637 

.901 

.995 

.991 

.996 

.05  (20) 

.662 

■  903 

.982 

.988 

1 .000 

.10  (10) 

.651 

CD 

VO 

CD 

.975 

1  .000 

1 .000 

.25  (6) 

.686 

.910 

1 .000 

1.000 

1 .000 

•50  (2) 

.750 

1 .000 

1 .000 

1 .000 

1  .000 

Properties  of  Truncated  Poisson 

Expected 

.632 

.896 

.999 

1 

1 .000 

1.000 

Ms-r  %  error  36.8  at  0.5  I0.lt  at  0.5 
Max  for  PK  £  0.5  36.8  at  0.5  10. <t  at  0.5 


0.1  at  0.5 
0. I  at  0.5 


Min  *  error  0.2  at  0.005  0.1  at  0.005  -2.5  at  0.10  -1.2  at  0.05  -0.6  at  0.025 

Min  for  PK  S  0.5  0.2  at  0.005  0.1  at  0.005  *2.5  at  0.10  -1.2  at  0.05  -0.6  at  0.025 

Avg  %  error  8.2  2.0  -0.6  -0.3  -0.2 

PK  S  0.5 


TABLE  A-12.  Comparison  of  Binomials  and  Truncated  Poisson  for  H  =  2. 


N 

1 

2 

5 

10 

20 

Expected  Kills 

i  for  the  Weighted  Binomials 

SLS/MSE 

.05  (60) 

.87 

1.67 

1.97 

1.96 

1  .98 

.10  (20) 

.88 

1.68 

1.96 

1.95 

2.00 

•25  (8) 

•  90 

1.69 

1.92 

2.00 

2.00 

■50  (6) 

.96 

1.56 

2.00 

2.00 

2.00 

.67  (3) 

.96 

1.65 

2.00 

2.00 

2.00 

1.0  (2) 

1 .00 

2.00 

2.00 

2.00 

2.00 

Properties  of  Truncated  Poisson 

1 

Expected 

.87 

1.66 

1.58 

2.00 

2.00 

Max  I  error  13-0  at  1.00  27.0  at  1.00 
Max  for  P|<  <  0.5  7-6  at  0.5  6.6  at  0.5 


1.0  at  1.00 
0  at  0.5 


Min  error 
Min  for  Pk  $  0.5 

Avg  %  error 
PK  S  0.5 


0  at  0.05 
0  at  0.05 


0.7  at  0.05  -3.1  at  0.25  -2.5  at  0.10  -1.0  at  0.05 
0.7  at  0.05  -3.1  at  0.25  -2.5  at  0.10  -1.0  at  0.05 
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TABLE  A-13.  Comparison  of  Binomials  and  Truncated  Poisson  for  H  =  5. 


N 

1 

2 

5 

10 

20 

Expected  Kills  for  the  Weighted  Binomials 

SLS/MSE 

.05  (100) 

•  99 

1.96 

6.14 

4.95 

4.88 

.10  (50) 

1  .00 

1.97 

4.14 

4.83 

4.82 

.25  (20) 

1  .00 

1.97 

4.04 

4.68 

5.00 

■50  (10) 

1  .00 

1 .98 

4.06 

5.00 

5.00 

•  67  (7.5) 

1  .00 

1.99 

4.23 

5.00 

5.00 

.83  (6) 

1 .00 

2.00 

4.50 

5.00 

5.00 

1.00  (5) 

1.00 

2.00 

5.00 

5.00 

5.00 

Properties  of  Truncated  Poisson 


Expected 

.99 

1.95 

4.12 

4.98 

5.00 

number 

kl  1  led 

Max  1  error 

1 .0  at  1 .00 

2.5  at  1.00 

17-6  at  1 .00 

0.4  at  1.00 

0 

Max  for  P|^  <  0.5 

1 .0  at  0.50 

1.5  at  0.5 

0.5  at  0.10 

0.4  at  0,5 

0 

Min  X  e  r  ro  r 

0  at  0.05 

0.5  at  0.05 

-2.0  at  0.25 

-6.4  at  0.25 

-3-7  at  0. 10 

Min  fo'  P|<  <  0. 5 

0  at  0.05 

0.5  at  0.05 

-2.0  at  0.25 

6.4  at  0.25 

-3-7  at  0. 10 

Avg  %  error 

0.8 

1 .0 

-0.6 

-2.4 

-1.5 

PK  <  0.5 

TABLE  A-14.  Comparison  of  Binomials  and  Truncated  Poisson  for  H  =  10. 


" — ..  N 

1 

2 

5 

10 

20 

Expected  Kills  for  the  Weighted  Binomials 

sls/mse 

.05  (200) 

1 .00 

2.00 

4.96 

8.78 

9-94 

.10  (100) 

1  .00 

2.00 

4.97 

8.74 

9.68 

■25  (40) 

1 .00 

2.00 

4.95 

8.40 

9.38 

•50  (20) 

1  .00 

2.00 

4.92 

8.31 

10.00 

.67  (15) 

1  .00 

2.00 

4.94 

8.57 

10.00 

.83  (12) 

1.00 

2.00 

4.97 

9.20 

10.00 

1.00  (10) 

1 .00 

2.00 

5.00 

10.00 

10.00 

Properties  of  Truncated  Poisson 


Expected 

number 

killed 

1.00 

2.00 

4.96 

8.75 

10.00 

Max  1  error 

0 

0 

0.8  at  1 .0 

12.5  at  1.00 

0 

Max  for  P«  <  0.5 

0 

0 

0.2  at  0.10 

0.3  at  0.05 

0 

Min  1  error 

0 

0 

-0.8  at  0.50 

-5-3  at  0.5 

-6.6  at  0.25 

Min  for  P«  S  0. 5 

0 

0 

-0.8  at  0.50 

-5-3  at  0.5 

-6.6  at  0.25 

Avg  X  error 

PK  S  0.5 

0 

0 

-0.2 

-2.3 

-2.6 
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